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ABSTRACT 

We performed spectroscopic X-ray observations of the eastern and northern 
regions of the Cygnus Loop with the ASCA observatory. The X-ray surface 
brightness of these regions shows a complex structure in the ROSAT all- sky sur- 
vey image. We carried out a spatially-resolved analysis for both regions and 
found that kT e did not increase toward the center region, but showed inhomo- 
geneous structures. Such variation cannot be explained by a blast wave model 
propagating into a homogeneous interstellar medium. We thus investigated the 
interaction between a blast wave and an interstellar cloud. Two major emis- 
sion mechanisms are plausible: a cloud evaporation model and a reflection shock 
model. In both regions, only a reflection shock model qualitatively explains our 
results. Our results suggest the existence of a large-scale interstellar cloud. We 
suppose that such a large-scale structure would be produced by a precursor. 

Subject headings: ISM: individual (Cygnus Loop) - supernova remnants — X- 
rays: Spectra 



1. INTRODUCTION 

The Cygnus Loop is the prototype of evolved shell-like supernova remnants (SNR). Its 
large apparent size and high surface brightness make it an ideal object for the study of the 
interaction between the supernova blast wave and the interstellar matter or the interstellar 
cloud in detail. Theoretical calculations of these interactions have been carried out by many 
authors (Sgro 1975; McKee & Cowie 1975; McKee & Ostriker 1977; Spitzer 1982; McKee 
1987; Klein, McKee & Colella 1994). 
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Various optical and ultraviolet observations of the Cygnus Loop have been performed to 
study the shock waves and their interactions. Cox (1972) suggested that the Cygnus Loop 
is in the radiative cooling phase of its evolution. X-ray observations, however, revealed that 
the Cygnus Loop is in the adiabatic phase and is still evolving (e.g. Ku et al.1984; Charles 
et al. 1985; Miyata et al. 1994). McKee & Cowie (1975) modeled the interaction between 
a blast wave and an interstellar cloud and were able to explain the discrepancies between 
the optically measured velocities of the filaments and the blast wave velocity inferred from 
X-ray observations. 

There are two models to explain the apparent optical and X-ray emission; an evaporation 
model and a reflection shock model. The evaporation model was proposed by McKee & Cowie 
(1975). They modeled a blast wave propagating into an inhomogeneous interstellar medium 
where there are many small clouds (cloudlets). In this model, the cloudlets are not directly 
heated by the shock wave but are evaporated due to thermal conduction by the ambient hot 
gas. This model was applied to both optical data (Fesen, Blair, & Kirshner 1982) and X-ray 
data (Tsunemi & Inoue 1980; Ku et al. 1984). It was, however, ruled out for the eastern 
region of the Cygnus Loop (e.g. Hester & Cox 1986, hereafter HC). Recent observations of 
[Ne V] emission at the eastern region also suggests that thermal evaporation does not play 
an important role (Szentgyorgyi et al. 2000). 

The reflection shock is important when the cloud size is large (~ lpc). Recent optical 
and X-ray data have been explained in the context of this model at various portions of the 
Cygnus Loop; the southeastern blob (SE blob) by Graham et al. 1995, the eastern region 
by HC and Raymond et al. 1988, and the northeastern (NE) region by Hester, Raymond & 
Blair 1994. 

Such a large scale cloud reminds us of a pre-existing cavity wall produced by the pre- 
cursor, as was first proposed by McCray & Snow (1979). Shull & Hippelein (1991) measured 
proper motions at 39 locations within the Cygnus Loop and found an asymmetric expan- 
sion of the shell. They also concluded the existence of a pre-existing cavity wall around 
the Cygnus Loop. This model is also supported by observations of coronal iron line emis- 
sion (Teske 1990), observations with IRAS (Braun & Strom 1986) and observations with 
the Rosat HRI (Levenson et al. 1997). Miyata & Tsunemi (1999, hereafter MT) found a 
large density discontinuity at the NE region. Combined with a large estimated mass of the 
progenitor star, they also supposed the existence of a pre-existing cavity. 

In this paper, we carried out spatially-resolved X-ray spectroscopy of the eastern and 
northern regions of the Cygnus Loop. Data analyzed here were obtained with the ASCA 
Observatory and the ROSAT Observatory. 
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2. OBSERVATIONS AND DATA CORRECTIONS 

We observed the eastern and northern regions of the Cygnus Loop in 1993-1994 with the 
ASCA observatory (Tanaka, Inoue, & Holt 1984). Figure 1 shows the location of SIS FOVs 
with black squares superimposed on an X-ray surface brightness map of the Cygnus Loop 
(Aschenbach & Leahy 1994). A detailed observational log is given in Table 1. We used the 
Faint-mode data obtained with the SIS to achieve the best spectral resolution (Yamashita 
et al. 1997). We excluded the data taken at an elevation angle below 5° from the night 
earth rim and 50° from the day earth rim, a geomagnetic cutoff rigidity lower than 6 GeV 
c -1 , and the region of the South Atlantic Anomaly. We also eliminated times when the 
pulse heights of corner pixels rapidly change. We then removed the hot and flickering pixels 
and corrected CTI, DFE, and Echo effects (Dotani et al. 1995; Dotani et al. 1997) in our 
data sets by using sispi and faint commands. Since we performed the observation of the 
northern region in the ASCA PV phase, the residual dark distribution (RDD) effect is not 
serious even for 4CCD mode data. However, the data of the eastern region observed in 1994 
were degraded due to the RDD effect (see detail in Dotani 1995). Therefore, we corrected 
the RDD effect for the eastern region using the frame data obtained near to the epoch of our 
observation. We retrieved the frame data taken for the calibration observation of Cyg X-l 
and anonymous sky (their sequence numbers are 18941107, 10021040 and 10021060) since 
they were observed in 1994 November. 

X-ray spectra accumulated in four regions are shown in Figure 2. We can confirm that 
the X-ray emission is dominated by thermal emission and emission lines from Mg and Si are 
clearly detected whereas emission lines from O, Ne, and Fe are weak. 

3. DATA ANALYSIS 

We investigated the spatial distribution of physical parameters using the SIS data in 
detail, using the same method as in MT in the northern region. In the northern region, we 
divided our FOV into 3' squares each of which corresponds to the half-power diameter of the 
XRT (Serlemitsos et al. 1995). Furthermore, we arranged the small squares such that each 
square overlapped adjacent squares by 2' resulting in 400 spectra. 

In the eastern region, however, the surface brightness is lower than that of the northern 
region as shown in Table 1. We divided our FOV into 5' in order to obtain statistics similar to 
those for the northern region. Each small region is arranged such that each square overlapped 
adjacent squares by 2' resulting in 49 spectra from each eastern region. 

Because the optical axes of SIS0 and SIS1 were different from each other, we employed 
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the ascatool library to adjust the optical axes. After adding the spectra of SIS-0 and 
SIS-1, we applied the single component non-equilibrium ionization (NEI) model coded by 
K. Masai (Masai 1984; Masai 1994) to each spectrum. Free parameters of the fitting were 
kT e , ionization parameter (log(r)), abundances of O, Ne, Mg, Si, Fe, and Ni, and emission 
measure (EM in cm~ 6 pc). Values for the reduced x 2 ranged from 1 to 2. 

3.1. Plasma Structure of the Eastern Region 

Figure 3 shows the spatial variation of kT c and log(r) for the eastern region. We find a 
complicated kT e structure revealed by the SIS. The values of kT c do not smoothly increase 
toward the center of the Cygnus Loop, but have local maxima at the crossing point of 
three observations (~ 0.6 keV) and the northern region of pos-EU (~ 0.8 keV). This kT e 
distribution cannot be explained with a blast wave model propagating into a homogeneous 
medium. On the contrary, the spatial variation of log(r) is relatively uniform with 10.3 ±0.4. 
We can confirm the hot plasma is still in the ionizing phase. 

3.1.1. Comparison with the ROSAT Image 

We retrieved the ROSAT PSPC data from the HEASARC /GSFC. The sequence number 
was RP500090. The X-ray surface brightness map obtained with the ROSAT PSPC is shown 
in Figure 4 (upper). Black squares show the FOVs of the SISs. Vignetting corrections were 
performed using the EXSAS/MIDAS software package. We see the instrumental structure 
due to the ribs of the window support of the PSPC which are seen only in a direction parallel 
to the direction of the telescope wobble. We found significant inhomogeneities in our FOVs. 
They could not be explained with a simple blast wave model. There is a bright clump (cl) 
at the PSPC center (a = 20 h 57 m 14 s ,5 = 31°1'55") with a radius of 2/7 (~ 0.6 pc) where 
we guess that the blast wave is moving into a dense region and has stagnated. The region 
named cl is the bright part of "XA" named in HC. The regions to the north and south of 
cl seem to have lower density, which means that the blast wave goes ahead of cl. There is 
another bright clump to the south of cl which is the SE blob. To the west of cl, we see a 
bright arc structure which is the left portion of pos-ER. 

The lower panel of Figure 4 shows the surface brightness map obtained with the SIS. 
Comparing these two images, we notice that cl, marked with a black circle, is not bright in 
the SIS image. The bright regions in the ASCA energy band (0.5—10 keV) are the northeast 
region of pos-ER and the southwestern region of pos-EU. The effective energy range of the 



- 5 - 



PSPC is 0.1-2.4 keV while that of the SIS is 0.5-10 keV. There is a big difference in 
the sensitivity below the C-K edge. Therefore, possible interpretations of the discrepancies 
between the ROSAT and ASCA images are 1) the presence of a low kT e component at cl, 2) 
the presence of a high kT e component at ASCA bright region, and/or 3) an extremely low 
Nn around cl. 

3.1.2. Band Ratio Map of the PSPC 

We constructed a band ratio map using the PSPC data by dividing the 0.44—2.0 keV 
band (R3— R7) image by the 0.11—0.28 keV band (Rl— R2) image, shown in Figure 5. We 
clearly see a soft X-ray emitting region along the shock front in the FOV. The softest region 
is at cl in this figure, shown in white. We also found a soft X-ray emitting region around cl 
as shown by a circle (hereafter referred to as c2). Around c2, a hard X-ray emitting region 
extends and wraps around c2 where kT e determined with the SIS is also high as shown in 
Figure 3. 

3.1.3. Combined Spectral Analysis 

As shown in Figure 5, a spectral variation is seen in the PSPC data. We analyzed 
the PSPC spectra at various positions within the FOV. We used xselect to extract the 
X-ray spectra, and used the response function pspcb_gain2_256 .rmf for the spectral fits. 
The effective areas of the XRT were calculated by pcarf , which included vignetting effects. 
Figure 5 shows the five regions (cl, rl— r4) where spectra were extracted. The spectrum 
obtained for each region is shown in Figure 6. The background spectrum was extracted 
from a blank region in the same data set. We clearly see strong soft X-ray emission below 
the C K-edge in all regions. This soft excess could account for the difference in the surface 
brightness maps between the ROSAT and ASCA energy bands. 

We extracted the SIS spectra from the same region as the PSPC (cl, rl— r4) and simul- 
taneously fitted the PSPC and SIS spectra. We applied the NEI model for these spectra. All 
parameters except the normalization were linked between the SIS and PSPC spectra. Free 
parameters were kT e , log(r), abundances of O, Ne, Mg, Si, and Fe, Ah, normalizations for 
both the SIS spectra and the PSPC spectra. We could not obtain statistically acceptable 
fits with a single component NEI model. In all PSPC spectra, we found excess emission 
below the C K-edge. We confirmed that such excess emission could not be reproduced with 
extremely low A H such as 10 19 cm -2 . We therefore added another component with a dif- 
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ferent kT e to reproduce the low energy region of the spectra. The low kT e component was 
mainly detected in the PSPC energy band. Due to the insufficient energy resolving power of 
the PSPC, we could not determine the individual abundance of heavy elements and log(r) 
for the low kT c component. We therefore assumed cosmic abundances and the condition 
of ionization equilibrium for the low kT c component. We obtained acceptable fits for all 
regions. Best fit curves are shown in Figure 6 and parameters are summarized in Table 2 
and 3. 

Including the ROSAT PSPC (C-band), we found a very low kT c component. The 
ASCA SIS is insensitive to such a low kT e component. The kT c values obtained for high kT e 
components are similar in all regions and also similar to those obtained at the NE region 
(MT). The values of log(r) are also similar in all regions whereas they are slightly higher 
than those in the NE region. 

Since the kT e value obtained for low kT e component is extremely low, the major part 
of its emission is below the PSPC band. Therefore, the kT e value obtained might contain 
some uncertainties because we only observe a small fraction of the spectrum. Moreover, the 
emissivity of such \ow-kT e plasma is very large (typically 1CT 22 erg s _1 cm -3 for plasma with 
cosmic abundances; Gehrels & Williams 1993), resulting to quickly cool down to ~ 10 5 K. 
Other narrow band studies, however, reveal the existence of such low kT e plasma. Images 
of [Fe X] and [Fe XIV] lines showed the existence of the cooler X-ray emitting plasma in the 
vicinity of the eastern region of the Cygnus Loop (Teske 1990; Teske & Kirshner 1985; Ballet 
et al. 1984). The image of the [Ne V] line showed bright filaments in the vicinity of cl, rl, 
and r2 regions (Szentgyorgyi et al. 2000) whereas it did not cover other regions. As [Ne V] 
is a unique probe of plasma with (1 — 6) x 10 5 K, we confirm that the extremely low kT c 
component is surely present at cl and rl. 

We should note that the metal abundances are depleted at cl and r4 while they are 
comparable to the cosmic values at rl, r2 and r3. This reminds us the spatial structure 
detected at the NE region of the Cygnus Loop. In the NE region, abundances of heavy 
elements are also depleted and shows an increase toward the inner region (MT). At the 
shell, a radius of > Q.9R S in the NE region where R s is the shock radius, they are constant 
(0.15~0.3 times cosmic values), whereas, behind the shell, they are 20—30% larger than 
those obtained at the shell. Therefore, abundances at cl and r4 are consistent with those 
obtained at the shell whereas those obtained at rl, r2, and r3 are qualitatively consistent 
with those obtained behind the shell with considering the statistical errors. 
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3.2. Plasma Structure of the Northern Region 

Figure 7 shows the spatial variations of kT c and log(r) in the pos-NN. As shown in 
Figure 1, a V-shaped structure is clearly seen over the SIS FOV. Values of kT c at the 
V-shaped structure (~ 0.3 keV) are lower than those in the neighboring region. This is 
consistent with results obtained with the ROSAT PSPC (Decourchelle et al. 1997). There 
are three high kT c regions marked in Figure 7. Such inhomogeneities in the kT e distribution 
cannot be explained by a simple blast wave model. The spatial variation of log(r) shows an 
anti-correlation with that of kT e . 

3.3. Comparison with the FPCS Results 

One of the most important purposes of observing the northern region is to compare 
the results with those of the FPCS (Vedder et al. 1986). The FPCS onboard the Einstein 
Observatory was the only available data with better energy resolution than ours and had a 
capability to resolve the triplet of O VII; forbidden (561 eV), inter-combination (569 eV), 
and resonance (574 eV) lines. 

To compare our results with those of the FPCS, we investigated the intensity of indi- 
vidual emission lines. We fitted the spectrum, extracted from the same region as Vedder et 
al. (1986), with a thermal bremsstrahlung continuum and then added a Gaussian line profile 
one by one until an acceptable fit was obtained (combination model as mentioned in Miyata 
et al 1998). Figure 8 shows the best fit curve and the residuals between the data and the 
model. An acceptable fit was obtained with a model that included 10 Gaussian line profiles. 
The best fit kT e and line intensity ratio of O VIII to O VII values are 0.31 ± 0.01 keV and 
0.91 ±0.06. 

Due to its insufficient energy resolution, the SIS cannot resolve the O VIII Ka line from 
the O VII K[3 line. We then calculated the line intensity ratio of the sum of O VII K(3 and 
O VIII to O VII based on the Table 1 in Vedder et al. (1986). Thus, the values obtained with 
the FPCS are kT c > 0.3 keV and a line ratio of 1.5 ± 0.6, respectively. The SIS results fully 
agree with the FPCS results to within the statistical uncertainties. We therefore confirm the 
consistency between the SIS and the FPCS. 
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4. INTERACTION BETWEEN THE BLAST WAVE AND INTERSTELLAR 

CLOUD 

As mentioned in section 3, a simple blast wave model is able to explain the inhomo- 
geneities neither in the eastern nor in the northern regions. We need to introduce interstellar 
clouds into the model in order to account for the inhomogeneities. The interaction of the 
shock wave with the interstellar clouds results in three kinds of shock; a shock traveling 
around the cloud, a shock transmitted into the cloud (hereafter referred to as the trans- 
mission shock), and a shock reflected from the cloud (hereafter referred to as the reflection 
shock). The nature of the interaction depends on the size and the density of the cloud. 



After the interaction, the reflection shock will propagate toward a direction opposite 
to the incident shock when the postincident shock flow is supersonic relative to the cloud. 
We assume radiation can be ignored and that the cloud is spherical. If the density contrast 
between the cloud and the ambient matter is not high enough to be treated as incompressive, 
some amount of energy will be released by the transmission shock, propagating into the cloud. 

The physical properties of the transmission shock as well as the reflection shock depend 
on the density contrast, a = n 4 /n , where is the density of the cloud. Sgro (1975) studied 
the interaction by considering the density contrast in the assumption of the strong shock 
limit for the blast wave (Mach number (M ) > 6). He solved the fluid equations for the 
mass, momentum, and energy conservation across the discontinuity between the reflection 
shock and the transmission shock. The geometry of the interaction is shown in Figure 9. We 
denote by subscripts 0, 1, 2, 3, and 4 the parameters of preshocked gas, postincident shocked 
gas, postreflected shocked gas, posttransmitted shocked gas, and incident cloud, respectively, 
as shown in this figure. Across the discontinuity, the pressure is constant, p 2 = Pz- When 
the postincident shock wave strikes the cloud, it would experience an overpressure relative 
to the ambient postincident shocked gas by a factor of (3. 



He found the following equations for temperature and velocity, with respect to the contact 
discontinuity, of both the reflection and transmission shock region as a function of postinci- 
dent shocked parameters. 



4.1. 



Reflection Shock Model 



n 3 vl = ftnivl 



(1) 



T 3 = r - 




(2) 



a 
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T 2 = (3) 

a r 

V3 = \ - Vi (4) 

V a 

* - {i-mfh • 

(6) 

where a r is the density contrast between the postincident shocked gas and the postreflected 
shocked gas, n 2 /ni, and can be related to a through 

3a r (4a r - 1) 

a ~ [{3a r (4-a r )} 1/2 - VE(a r -l)]* 
Figure 10 shows the characteristics of the postreflected shocked region as a function of a. 
This figure is consistent with Figure 17 of Hester et al. (1994) even though their approach 
of calculation differs from ours. 



4.2. Cloud Evaporation Model 

When the postincident shock wave strikes the cloud, it would experience an overpressure 
relative to the ambient postincident shocked gas by a factor of (3. Values of temperature 
and velocity for the transmission shock can be determined by equations (1,2). The time 
required for the transmission shock to travel through the cloud is 10 6 i?4/t>3 yr, where R4 is 
the radius of the cloud in pc and V3 is the velocity of the transmission shock in km s -1 . By 
equating the cooling time of the posttransmitted shocked gas to the time required for the 
transmission shock to travel through the cloud, a critical cloud density, n c , is defined. Sgro 
(1975) calculated n c and the cooling time, r c , of the posttransmitted shocked gas due to the 
radiation as 

r c = 4.6 x 10 6 (^) % 2 vt yr, (9) 

where e is the emissivity of thin thermal gas, which is approximately 1.33x 10 _16 T _1 ergs 
cm" 3 s" 1 in the range of 10 5 < T < 10 7 K. 

For clouds with density greater than n c , the posttransmitted shocked gas cools rapidly 
and clouds begin to evaporate into a hot tenuous medium. Note that (3 depends on the 
density contrast, a = n c /n , so that we should pay attention to calculate n c . 
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Figure 11 shows the maximum cloud radius which can be heated up by the blast wave. 
f3 is calculated for each density contrast and Ti is fixed to be 0.21 keV. This graph allows us 
to determine n c for a given n and R c . 

After the cloud is surrounded by the tenuous medium, it will be gradually heated up 
due to thermal conduction. The thermal conductivity of a fully ionized Hydrogen plasma is 
given by (Spitzer 1962) 

« = 1-94 x 10 11 (y^) 5/2 (^) 1 ergs s^K-W 1 , (10) 

where In A is the Coulomb logarithm which we set to lnA=30 since it is insensitive to T c 
and density. We set the differences in temperature between the cloud and the surrounding 
medium to be AT, so that the heat flux is 

q = -nAT . (11) 

The heat flux is proportional to X/Lt, where A is the mean free path for electron energy 
exchange and L T is the temperature scale height T/\AT\. X is found from 

, 3kT^ 1/2 

A = t, 
= 0.2 i 



m 



IwO ho" pc • (12) 



10 7 Ky Vlcnr 

The classical model of thermal conductivity assumes that the mean free path, A, is 
shorter than Lt, which is a typical cloud size. However, for many situations of interest to 
our model, the mean free path of electrons is not shorter than the radius of the evaporating 
cloud. When A becomes comparable to or greater than Lt, the heat flux is no longer 
equal to — kAT. This effect is referred to as saturation. Cowie & McKee (1977) studied 
'saturated' thermal conduction in detail. As to the spherical cloud, the timescale for complete 
evaporation can be written by (Cowie & McKee 1977) 

, . . . R 4 \ 2 / Ti V 5/2 /In A 
t. ' 5 ' 



u evp 



5. SHOCK STRUCTURE OF THE EASTERN REGION 
5.1. Density of Bright Clump 



As shown in Figure 5, we found two soft X-ray bright regions in the eastern region; cl 
and c2. Inside c2, we detected low kT e components at every position (cl and rl— r4). Based 
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on these results, we calculated the density of these regions. We assumed a homogeneous 
spherical cloud for c2 with a radius of 12' m 2.7 pc for the low kT c components. The 
calculated volumes and densities are summarized in Table 4. If we assumed pressure balance 
between the low kT e and high kT c components, we could calculate the density and volume 
for the high kT e components as shown in Table 4. 

The obtained density for the low kT e component is roughly 30 cm -3 at every portion of 
the eastern region. The density for the high kT e component is roughly 5 times higher than 
that obtained at the NE region as found in MT. 

We also calculated the thermal pressure at each region shown in Table 4. HC estimated 
the thermal pressure at the "XA"to be (1.1 — 3.8) x 10 7 cm -3 K, which is consistent with 
our results. 



5.2. Emission Mechanism of the Eastern Region 

Because a simple blast wave model can not explain the inhomogeneities of the eastern 
region, we investigate emission mechanism models which can hold for the overall X-ray 
features of the eastern region. 

5.2.1. Locally Large Preshock Density 

We first consider an inhomogeneous ISM. Specifically, we consider an enhancement of 
the preshock density. To distinguish it from the interstellar cloud, we assume that the 
enhanced density is lower than n c . In this case, such an enhancement will be smoothed out 
by the blast wave, resulting in an enhancement of the X-ray emission. This could account 
for the bright X-ray appearance at the eastern region. 

Based on this model, kT e of the bright region is expected to be lower than those of 
neighboring regions. The high kT e region can be considered as the postincident shocked 
region. This means that we can apply a simple blast wave model to the high kT e region as 
a first order approximation. Based on the Sedov model (Sedov 1959), kT c just behind the 
blast wave is proportional to n™ 2 / 5 . For the high kT e components, n e is roughly 6 cm -3 . 
Assuming the strong shock limit, no is calculated to be 1.5 cm" 3 . This value is roughly 6 
times higher than that of the NE region (MT). Considering the kT e value at the NE region, 
kT e at the eastern region should be (1.5/0.25)" 2 / 5 x 0.3 0.15 keV. This value is roughly 
two times lower than that obtained at c2 region as shown in Table 4. Therefore, the kT c 
value obtained for the high kT c component is too high to be explained using the simple blast 
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wave model. 

Thus, it seems difficult to explain the current bright X-ray features as resulting from a 
shocked region in which the density is higher than that of the neighboring region. 

5.2.2. Cloud Evaporation Model 

We found region c2 to be bright in soft X-rays. As shown in Table 4, the density of 
this region is quite high compared with other regions (0.25 cm" 3 at the NE region). If we 
considered c2 as an interstellar cloud, cloud evaporation might play an important role in 
its high density. Based on the cloud evaporation model, however, high kT e components are 
considered as the emission from the postincident shocked region just as the same case of 
the previous model. Thus, we met with a similar difficulty in trying to explain the X-ray 
properties by the cloud evaporation model. 

We next calculated the cooling time of the posttransmitted shocked gas. For clouds 
with density greater than n c , the posttransmitted shocked gas rapidly cools down to 10 5 
K, and there is no X-ray emission. Based on the equation (8), we calculated the maximum 
cloud thickness, R4, to be heated by the transmission shock to emit X-rays. Figure 11 shows 
the maximum cloud radius plotted as a function of a and no- Vertical dashed lines show 
the permitted area for a (inferred in the subsequent subsection) while horizontal solid lines 
show the permitted area for hq. This figure shows that the transmission shock wave will cool 
down to below the X-ray emitting temperature when it moves into a cloud roughly 0.02—3 
pc in size. This indicates that the transmission shock can heat up roughly half of c2 for 
which the radius is roughly 2.7 pc and it cannot pass through c2 since it will cool down to 
below the X-ray emitting temperature. Therefore, the cloud evaporation model can be ruled 
out as an explanation for the hard X-ray emitting gas, although it might account for the 
\ow-kT e component for future observations. 

5.2.3. Reflection and Transmission Shock Model 

As described in the previous section, three kinds of shock waves might occur when 
the blast wave interacts with the cloud. We suppose that the posttransmitted and the 
postreflected shocked regions correspond to the soft and hard X-ray regions, respectively, in 
case of the eastern region. The reflection shock propagates into the postincident shocked 
region and provides an additional heating. The transmission shock propagates into the high 
density cloud and heats it up to X-ray emitting temperatures if the density contrast between 
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the cloud and the ambient matter is not too high, accounting for the bright X-rays seen in 
the ROSAT band. 

Based on this model, the ratio of kT c for the postreflected shocked region to that for 
the posttransmitted shocked region depends on the density contrast. Based on Figure 10 
we calculated the ratio for kT c and n e between the postreflected shocked region and the 
posttransmitted shocked region, shown in Figure 12. We also plotted our results in Figure 12 
for cl and rl— r4. The density contrast, a, can be determined based both on the kT e ratios 
and the n e ratios as shown in Table 5. We can obtain a consistent value of a ~ 10 for all 
regions even though the statistical uncertainties are large for the n e ratios. Therefore, based 
on the reflection shock model, we can quantitatively explain the observed values of kT e for 
high and low kT c components with an interstellar cloud for which the density is roughly an 
order of magnitude larger than that of the ambient medium. 

We can calculate the various physical parameters of the postreflected and posttransmit- 
ted shocked region based on Figure 12. Table 5 shows values of n 2 /n 1 (= a r ) and T 2 /Ti 
for each region. Based on these results, n e and kT c at the postreflected shocked region is 
roughly 1.8 and 1.5 times larger than those at the postincident shocked region. Therefore, 
n e and kT e at the postincident shocked region are ~ 2.8 cur 3 and ~ 0.21 keV. This implies 
that no at the eastern region is 0.7 cm -3 . This value is roughly three times larger than the 
mean value (Ku et al. 1984). This might be a region of enhanced density, a 'tail' of the 
cloud, extending around c2 as suggested by HC. 

If we compare this density with that obtained at the NE region, kT e at the eastern 
region is found to be (0.7/0.25) -2 / 5 x 0.3 « 0.20 keV. This value is in good agreement with 
our result. 

Thus, we conclude that the reflection shock and transmission shock can explain the 
X-ray properties observed both with ROSAT and with ASCA. 

It is interesting to note that the density contrast obtained at the eastern region is 
consistent with that obtained at the NE region (MT). Combining these results, we suppose 
the existence of a large-scale high density region from the NE toward the eastern region of 
the Cygnus Loop. 
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6. SHOCK STRUCTURE OF THE NORTHERN REGION 

6.1. Complicated X-ray Surface Brightness and kT c Map 

The X-ray image of the pos-NN shows a peculiar V-shaped feature, as shown in Figure 1. 
This V-shaped region is the brightest among the regions with the same radius from the center. 
The kT e distribution of the pos-NN is not uniform as shown in Figure 7. We found that 
three regions showed kT c (0.5—0.6 keV) higher than other regions (0.3—0.4 keV). We refer to 
these three regions as hi— h3 as marked in Figure 7 with green lines. Comparing the surface 
brightness map with kT c map, we noticed a strong anti-correlation between them. 

6.2. Emission Mechanism of the Northern Region 

We centered our FOV on the brightest part in the northern region which was about 20' 
inside the shock front. As mentioned in section 5.2, we expected that the brightest region 
was roughly 4' inside the shock front if we assumed the blast wave to be propagating into a 
homogeneous medium. The apparent structure of the pos-NN is unlikely to be explained by 
such a simple blast wave model. 

As mentioned in the previous section, the pos-NN is the brightest among the regions 
with the same radius from the center. One of the following two scenarios may explain the 
emission mechanism; the high density medium may be like an interstellar cloud (evaporation 
or reflection shock heating) or the shock front structure may be due to a projection effect. 
In the subsequent section, we will investigate the possible explanation for the current X-ray 
features observed both with ROSAT and with ASCA. 



6.3. Shock Front Structures due to the Projection Effect 

The mean kT e around the V-shaped region is roughly 0.4 keV. It is similar to that 
obtained at the NE region (MT). However, the surface brightness is much higher than regions 
with the same radius in the direction of the NE region. If the V-shaped region comes from 
the shock front due to a projection effect, the mean density of the pos-NN must be higher 
assuming the pressure equilibrium. Therefore, we can expect a lower kT c than that of the 
NE region assuming the pressure equilibrium. However, we found kT e to be similar to each 
other. We conclude that the bright X-ray appearance does not come from a projection effect 
of the shock front structures. 
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6.4. Cloud Evaporation Model 

We next investigate a cloud evaporation model. After the interstellar cloud is engulfed 
by the shock wave, it begins to evaporate when the cloud density is higher than the critical 
density for direct shock heating. The value of kT c of the shocked cloud (kT 3 ) depends on the 
density contrast between the cloud and the ISM. The value of kT 3 must be lower than that 
of the postincident shocked region. Assuming that the cloud is sitting at the center of our 
FOV, we can expect bright X-ray emission from the cloud. This picture resembles our results 
at the point of 1) the bright and low kT e emission at the center of our FOV and 2) high kT e 
emission at h2 (and possibly h3). Generally speaking, the value of kT e decreases toward the 
shock front where the density is the highest. Therefore, hi is considered as the postincident 
shocked region shocked by the blast wave passing over the cloud. We then investigated 
whether or not we could explain the X-ray feature of hi by a model that included direct 
shock heating. 

The angular distance from the center of the Cygnus Loop to hi is roughly 72'. The 
value of kT e obtained at the region with the same radius in the direction of the NE region 
is roughly 0.4 keV (MT). If such a high kT e medium at hi is the shock heated ISM, we can 
estimate the density by using the Sedov model. Since kT e of the shock wave is proportional 
to n~ 2/5 , n of the pos-NN is expected to be (0.6/0.4)" 5 / 2 x 0.25 « 0.09 cm" 3 . The surface 
brightness of hi is lower by a factor of ~ 3 than that of the NE region, which is consistent 
with the low density derived above. However, based on a similarity solution, the shock radius 
is proportional to n -1 / 5 . The current location of the shock front passing over the cloud is 
expected to be (0.25/0.07) 1/5 « 1.3 farther away from the center than that at the NE region. 
However, it is not the case as shown in Figure 1. Therefore, it seems difficult to explain the 
current X-ray feature of hi with the model of shock heating by the incident shock wave. 

Therefore, even if we take into account not only a single cloud but also multiple clouds, 
we expect that the cloud evaporation model would be ruled out to explain the pos-NN. 

6.5. Reflection Shock Model 

We found that neither the direct shock heating model nor the cloud evaporation model 
could account for the pos-NN of the Cygnus Loop. We next consider the reflection shock 
heating model. 

There are several configurations for the location of the cloud. Among them, we choose 
two favorable configurations as discussed in the subsequent sections. 
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6.5.1. Cloud at the Northern Region in Our FOV (hi) 

We first examine the case that a cloud exists at hi. We can interpret the structure of 
the V-shaped region as the postreflected shocked region as shown in Figure 9. Regions of 
h2 and h3 are expected to be the postincident shocked region. In this case, we expect kT c 
at the V-shaped region to be higher than those of h2 and h3. However, this is not the case. 
Therefore, this is unlikely to be explained with the reflection shock model. 

The value of kT c of hi is much higher than those at other regions at the same angular 
distance from the center. Such high kT e is expected in the postreflected shocked region. If a 
cloud is at hi, we will expect the emission from the cloud itself at hi either in at soft X-ray 
or optical wavelengths. 

The value of kT c of the posttransmitted shocked region depends on the density contrast. 
For a small density contrast, we expect that a low kT e component is bright in the soft X-ray 
energy band (ROSAT band). On the contrary, for the case of a high density contrast, kT c 
of the posttransmitted shocked region decreases so that there is no X-ray emission, so that 
we then expect some optical counterparts to be seen at the position of the cloud. Based on 
the kT e map obtained with the ROSAT all-sky survey (Aschenbach & Leahy 1994), kT e at 
hi region is higher than those at its neighboring regions. This will deny the existence of the 
soft X-ray emitting component. 

Figure 13 shows the optical image of the pos-NN of the Cygnus Loop which is a digital 
sky survey image retrieved from the Skyview program, supported by HEASARC/GSFC. 
There are some optical counterparts seen. However, we find no emission around the hi 
region. Therefore, we conclude that there is not a low kT e component at regions neighboring 
hi either in the soft X-ray region or at optical wavelengths. 

Thus, this configuration cannot explain either the low kT e of the V-shaped region or the 
lack of low kT e component at hi. 

6.5.2. Cloud at the Center of Our FOV 

We next consider a cloud at the center of our FOV. After the interaction between 
the blast wave and the cloud, the blast wave passes over the cloud. The reflection shock 
propagates into the postincident shocked region while the transmission shock propagates 
into the cloud. The current geometry we suppose for the pos-NN is schematically shown in 
Figure 14. 

In this configuration, we directly observe the postreflected shocked region with a higher 
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kT e among the neighboring mediums at two different regions. Between these two regions, 
we observe the overlapping of a postreflected shocked region and a posttransmitted shocked 
region. Because the density of both regions is high, we expect substantially bright X-ray 
emission at the overlapping region. The value of kT e at the overlapping region is intermediate 
between the low kT e (kT 3 ) the high kT c (kT 2 ). 

If we apply this configuration to our data, we suppose that hi— h3 come from the 
postreflected shocked region as marked by upper and lower broken lines. Between hi and 
h2, the bright V-shaped region can be interpreted as the overlapping region. Therefore, we 
favor this picture since this picture can explain both the high X-ray surface brightness seen 
at the center of the pos-NN obtained with ROSAT and the kT e distribution determined with 
the ASCA SIS. 

We extracted spectra at hi— h3 regions and fitted them with a single component NEI 
model. However, we failed to fit them with a single component model. As shown in Figure 
7 there are still spectral variations in each region, possibly because there are some amount 
of contamination from the low-/cT e component and also the model we proposed is too much 
simplified. However, we believe the overall plasma structure can be explained with the 
current configuration. 

We estimate other physical parameters based on this model assuming a spherical cloud. 
The separation between hi and h2 roughly corresponds to the diameter of the cloud. The 
cloud is then at (a = 20 h 49 m 57 s ,5 = 32°3'33") with a radius of roughly 7' ~ 1.6 pc. The 
values obtained for kT c at hi— h3 are considered as T 2 whereas that of V-shaped region (~ 
0.4 keV) stands for the mean value of T 2 and T3. Therefore, we extracted the X-ray spectrum 
at the V-shaped region and fitted it with a two component NEI model with different kT e . 
The values obtained for kT c are ~ 0.2 and ~ 0.5 keV. The value of the high kT c component 
agrees reasonably with that obtained at hi— h3. Thus, T 3 can be considered as ~ 0.2 keV. 

The ratio of T 2 to T3 depends on the density contrast as discussed in section 4. Based 
on Figure 12, we obtain the density contrast a to be m 5. The cloud density is estimated to 
be roughly 1 cm' 3 assuming the ambient density of 0.2 cm -3 . For a = 5, we obtain kT c of 
postincident shocked region Ti of 0.45 keV and (5—2. 7\ obtained is similar to that of the 
region with the same angular distance from the center to the direction of the NE region. 

Figure 15 shows the maximum radius of the cloud which can be heated up directly by 
the shock wave, in the a — T\ plane. The vertical solid line represents the value of a at the 
pos-NN region. As shown in this figure, the cloud with size of < 1.6 pc will be destroyed by 
the blast wave with kT\ > 0.2 keV. On the contrary, the value of kT\ we obtained is 0.45 
keV. Such a cloud cannot survive during passage of the blast wave. Therefore, we can rule 
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out the cloud evaporation model to explain the bright X-ray features seen at pos-NN. 

In this way, the reflection shock model can well explain the observed X-ray features with 
a simple geometry as shown in Figure 14. 

7. CONCLUSION 

We have performed high quality X-ray observations of the Cygnus Loop using the ASCA 
Observatory and the ROSAT Observatory. We have found that the X-ray surface brightness 
and kT e significantly vary inside our FOV. They cannot be explained with a simple blast 
wave model. We carried out calculations for the interaction between the blast wave and the 
interstellar cloud in order to determine the characteristics of the postincident shocked gas, the 
posttransmitted shocked gas, and the reflection shocked gas. We applied the evaporation 
model as well as the reflection (and transmission) shock model for the eastern and the 
northern regions. We have difficulties explaining the X-ray nature with the evaporation 
model for both regions. Only the reflection shock model can quantitatively explain the 
observed kT e and the X-ray surface brightness. 

The cloud size is difficult to estimate. However, it is at least 1 pc in size for both regions. 
Such a large cloud is sitting in the NE region as well (MT). Combining all of these results, 
we conclude that a cloud of large scale is covering at least one third of the Cygnus Loop. 
This fact strongly supports the existence of a pre-existing cavity wall, possibly created by 
the precursor. 
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Table 1. OBSERVATIONAL LOG OF THE CYGNUS LOOP 



Date 


Sequence Number 


Coordinate (J2000) 


Position 


CCD mode" 


Exposure (ks) 


Intensity (c s 1 ) 


1993 Oct 18 . . . 


50008010 


20 h 50 m 2. s 4, 32 d 3' 4" 


pos-NN 


0123 / 2301 


5.5 / 6.5 


12.44 / 12.13 


1994 Nov 01 . . . 


52003020 


20 h 57 m 7. s 2, 31 d 14' 56" 


pos-EU 


0123 / 2301 


3.7 / 5.1 


6.60 / 5.78 


1994 Nov 01-02 


52003050 


20 h 57 m 22. s 3, 30 d 53' 13" 


pos-ED 


0123 / 2301 


5.9 / 6.5 


3.02 / 2.73 


1994 Nov 02 . . . 


52003000 


20 h 55 m 30. s 0, 30 d 56' 49" 


pos-ER 


0123 / 2301 


2.9 / 4.3 


5.16 / 4.00 



a CCD mode shows the CCD ID number used in the observation. 



Note. — The columns which have two values divided by '/' show the value of SIS-0 and SIS-1, respectively. 
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Table 2. Spectral Features at the Eastern Region with the Combined 
Spectral Analysis of the ASCA SIS and the ROSAT PSPC 



Region kT e -l log(r) EM-1 kT c -2 EM-2 AT H reduced \ 2 (dof) 

(keV) (cm- 6 pc) (eV) (cm- 6 pc) (10 20 cm" 2 ) 

el 0.27t°;o| 10.91q 4 20±f p 43±? (1.0±°;g)xl0 4 8.8 ± 0.7 1.2(116) 

rl 0.34±0.04 11.4 ±0.3 9±| 54.ll°^ (3.0±^)xl0 3 8.6±g-| 1.2(144) 

r2 0.34±°,'° 4 H.9lo3 3 ±2 67±f 3 (4^°')xl0 2 7.3±j'| 1.3(143) 

r3 0.3lt°''°| H.7t° : | 2± 6 54± 4 (3^°)xl0 3 8.7±^ 1.3(141) 

r4 0.32t°'' 66 10.8±J'| 12.31^ A7±\° (4±|)xl0 3 9±£' 1.1(121) 
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Table 3. Abundances of Heavy Elements Obtained at the Eastern Region 



Region C, N, O Ne Mg Si Fe 



cl 0.2l°i 0.31S1 0.25^3 <10 0.2t° i 

n o.8±8:| 1.71™ i.o±8:i o.9±8:g o.5oi° : g 

r2 o.s±S:| i.o^ o.9±»:S i-3^ o.5±S:i 

r3 6^ < 10 6t 2 5 < 6 2.8l°j 

r4 0.26^ 0.3±S;f < 0.3 <1 0.19±g;g 
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Table 4. Volume, Density, and Pressure of the Eastern Region 



Region 


High kT c 


Low kT c 


Pressure (p/k) 




Density 


Volume 


Density 


Volume 






(cm -3 ) 


(pc 3 ) 


(cm -3 ) 


(pc 3 ) 


(10 7 cm" 3 K) 


cl 




0.50 


40 ± 10 


7.4 


4±1 


rl 


4.6 ±0.6 


0.45 


29+ 3 


3.7 


o fi +0.3 


r2 
r3 
r4 


i-9tS:? 

°-3 


0.79 
0.073 
0.56 


w+l 1 

30l 30 

OQ + 30 

DU -20 


3.7 
3.7 
3.7 


o+3 
J -2 
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table 5. characteristics of the reflection and the transmission shocked 

Gas at the Eastern Region 



Region 


a (based on kT e ratio) 


a (based on n e ratio) 


n 2 /ni (a r ) a 


T 2 /T 1 a 


cl 
rl 
r2 
r3 
r4 


11 ±4 

10 ±5 
12±2 


11+? 

nil 5 

> 3 

> 4 

> 5 


1 8+ - 1 
1.8 ±0.1 

1.7 ±0.1 
1 7+ - 1 

L - ' -0.2 

1.8 ±0.1 


1.5 ±0.1 

i -°+-0.02 

1.5 ±0.1 

1.5 ±0.1 

1.6 ±0.1 


a These values are calculated using 


; a extracted from kT e 


ratio. 
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FIGURE CAPTIONS 

Fig. 1. — X-ray surface brightness map of the Cygnus Loop obtained during the ROSAT 
all-sky survey (Aschenbach & Leahy 1994). The black squares show the FOVs of the eastern 
and northern portion as observed with ASCA. 

Fig. 2. — X-ray spectra obtained at pos-NN (a), pos-EU (b), pos-ED (c), pos-ER (d). X-ray 
emission is dominated by thermal emission and emission lines from Mg and Si are clearly 
resolved. Line identifications are also shown in (a). 

Fig. 3. — (a) kT c and (b) log(r) maps of the eastern region obtained with the SIS. North is 
up and east is left. 

Fig. 4. — X-ray surface brightness map of the eastern region obtained with the ROSAT 
PSPC (upper) and the ASCA SIS (lower). Black squares show the FOVs of the SISs and a 
black circle marks the bright clump (cl). The difference between the surface brightness in 
the PSPC and that in the SIS energy bands is clearly seen which is due to the difference of 
the energy bands. 

Fig. 5. — Band ratio map of the eastern region of the Cygnus Loop. This map is extracted 
as (0.44-2.0 keV) / (0.11-0.28 keV). The region of cl and the soft X-ray emitting region 
(c2) are encircled by small and large circles, respectively. The regions where spectra (rl— r4) 
are extracted are shown. 

Fig. 6.— ASCA SIS and ROSAT PSPC spectra at cl (a) and rl-r4 (b-e) regions. Best fit 
curves are shown with solid lines. 

Fig. 7. — Same as Figure 3 but of the pos-NN. Three high kT e regions are marked by green 
circles, named hi— h3. 
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Fig. 8. — X-ray spectrum obtained from the portion of the remnant observed with the FPCS. 
Each component is also shown by a solid curve. 

Fig. 9. — Geometry of the reflection shock. Upper panel shows the cross section of the 
interaction with the cloud. Lower panel shows the density structure along the cross section 
of the cloud. 

Fig. 10. — Physical parameters for the reflection shock as a function of the density contrast 
between the cloud and the interstellar medium. 

Fig. 11. — Maximum cloud thickness which can be heated up by the blast wave is plotted 
versus a and n . The value of Ti is assumed to 0.21 keV. Contours are shown in a loga- 
rithmically scaled from 1 x 10~ 4 to 1 x 10 4 pc. The cloud size of c2 (2.7 pc) is also shown. 
Vertical dashed lines show the permitted area for a in the eastern region while horizontal 
solid lines show the permitted area for n . 

Fig. 12. — Ratios of the temperatures and of the densities between the postreflected and 
posttransmitted shocked regions. We should note that the density ratios plotted for r2, r3, 
and r4 only show the upper limits. 

Fig. 13. — Digital sky survey image of the pos-NN of the Cygnus Loop obtained with the 
Sky View. A black square shows the FOV of the SIS. 

Fig. 14. — Geometry of the interaction between the blast wave and the interstellar cloud 
which accounts for the X-ray appearance of the pos-NN. 

Fig. 15. — Same as Figure 11 but the vertical axis shows 7\. Ambient gas density is fixed 
to 0.2 cm~ 3 . Contours are shown in a logarithmically scaled and the cloud size of 1.6 pc is 
also shown. Vertical solid line represents the value of a obtained at the pos-NN. 
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